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Abstract 

We examine the Stillinger- Weber analysis of the potential energy landscape for its stability and conclude that it does not provide 
a stable description of the system as the free energy slope and curvature vanish simultaneously. An alternative analysis developed 
recently by us involving complexity provides a stable description with complexity a monotonic increasing function of temperature. 



It is well known that most supercooled liquids (SCL) be- 
come viscous when their configurational entropy becomes 
negligible as they are cooled, provided the correspond- 
ing crystal is not allowed to nucleate. Our current un- 
derstanding of glassy behavior is still far from complete, 
even after many decades of continuous investigation. In 
order to better understand the flow properties of viscous 
fluids, Goldstein proposed the potential energy landscape 
(PEL) picture using classical canonical ensemble (H 01 to 
qualitatively discuss an interesting but sufficiently tractable 
scheme to study SCL and the glassy states by drawing at- 
tention to the potential energy minima (having the energy 
£"), to be called basin minima (BM) in the following. The 
landscape picture with its BM's plays a pivotal role not 
only in the thermodynamics of viscous fluids at low tem- 
peratures but in many disparate fields like glasses, proteins 
and clusters |3], and has established itself as an important 
thermodynamic approach in theoretical physics. Thus, it is 
highly desirable to understand the significance of this ap- 
proach. Stillinger and Weber (SW) extended the work of 
Goldstein to higher temperatures by carrying out a formal 
analysis in terms of the minima energies £ |4, 5]. Their 
analysis has given rise to a considerable amount of litera- 
ture in recent years; for a partial list, see |0.|5l0,0|. Many 
of the numerous numerical evidence \1\ appear to be con- 
sistent with Goldstein's seminal ideas 1 1]. 

In this work, we study the stability of the SW analysis, 
which seems not to have been investigated in the litera- 
ture. There are two different conditions for the stability of 
a thermodynamic theory |8]. The first one is the vanish- 
ing of the slope of the free energy function and is com- 
monly discussed in the literature. The other condition is 
of a strictly positive curvature of the free energy function 
at the point where the first condition is met. This does not 
appear to have been ever discussed. To our surprise, we 
find that the free energy function in the SW analysis has a 
zero curvature. Thus, the SW analysis does not give rise to 
a stable description of the system and must be replaced by 
other self-consistent approaches. We have recently devel- 
oped such an approach fl fioll . which borrows the concept 
of complexity developed for spin glasses The new 

analysis has no problem with stability and is consistent. 

Conventional Approach. The canonical PF Z{T) for a 



system of N particles in a volume V is 



(1) 



where W{E) represents the number of configurations of 
energy E and defines the microcanonical entropy S{E) = 
In W{E), and /3 is the inverse temperature in the units of 
the Boltzmann constant. The value of Z(T) for a macro- 
scopic system, which is what we consider here, is de- 
termined by the dominant term in the sum, which is lo- 
cated at the equilibrium energy E = E[T): Z{T) = 
W[E] exp{—(3E). The determination of E for a macro- 
scopic system is simplified by noting that E is almost a 
continuous variable for a macroscopic system. In terms 
of S{E), E is given by the location of the minimum of 
the free energy function F[T, E) = E — TS{E) at fixed 
T. In equilibrium, the entropy S{T) = S{E = E) and 
free energy F{T) = F{T,E) = E{T) - TS{T) be- 
come functions only of T. The conditions for the minimum 
are [dF{T, E)/dE]T = 0, and [d^F{T, E)/dE\ > 
leading to 



[dS{E)/dE]^^ 



E=E(T) 



P, dS(T)/dT > 0, (2) 



which are always satisfied because of a non-negative heat 
capacity. With the use of (|2j, we immediately conclude 



T{dS{T)/dT) = dE{T)/dT, 



(3) 



which is consistent with the first law of thermodynamics at 
constant V, and A^. 

At a given temperature T, only those configurations that 
have the energy E = E {ox within a narrow width around 
it, depending on the heat capacity; we will neglect this 
width here) determine the thermodynamics through the en- 
tropy S{T). All energies other than E and, therefore, all 
configurations not included in W{E) are irrelevant at T. 
Thus, thermodynamics is highly selective. This will remain 
true even in the landscape picture, where the equilibrium 
states will have the same energy regardless of which basin 
they belong to. Thus, E = E(T) cannot depend explicitly 
on the basin minima energy £:{dE / d£)T = 0. 

Goldstein's Approximate Analysis. In his analysis, 
Goldstein has listed two conjectures that were common in 
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the field at the time: the basin PF Z\,{T) is (i) indepen- 
dent of the basin's minimum energy £, and (ii) insensitive 
to the basins being explored. Utilizing these assumptions, 
Goldstein has expressed the PF as a product ij^ of the basin 
and BM PF' 

Z(T) = Zb(T)ZBM(T); (4) 

here z\, for a given basin is defined by considering shifted 
energies E — £ with respect to the minimum energy £ of 
that basin; see also Schulz lfl3ll . Goldstein has emphasized 
that basin anharmonicity or the curvature at its minimum 
lfl4ll may be very important. These are included in Zb, so 
that it is determined by the entire basin topology. Accord- 
ing to Goldstein, all equilibrium basins have the same equi- 
librium basin free energy fh{T) = —Tin z^,. The BM-PF 
is defined Hill as 

Zbm{T) = Y,N^M{e)e-^' . (5) 

£ 

Here, A'^bm(^) represents the number of basins whose BM 
are at energy £. The equilibrium BM energy £ = £{T) is 
the value of £ at which the summand in (O is maximum. 
The conditions for the maximum in terms of the BM en- 
tropy Sbm{£) = IniVBM('f) are given by 

[dSBM{£)ld£],^j = d£/dT > 0, (6) 

which are the standard conditions of equilibrium; compare 
with Thus, the analysis is completely stable in this 
approximation. It is clear that the BM description pro- 
posed by Goldstein ensures that £ is a monotonic increas- 
ing function of T. Since this approximation is expected to 
be good at low temperatures, we expect £ to be monotonic 
increasing there. But it need not be true at all tempera- 
tures as shown below, thereby limiting the usefulness of 
the BM-description at all temperatures that has been for- 
mally adopted by Stillinger and Weber to which we now 
turn. 

SW Analysis. A basin is indexed by j, and the low- 
est and highest basin energies are denoted by £j , and £j , 
respectively, so that the basin does not exist outside the 
energy range AjE = i£j,£j). Let Wj{E) {E G /\jE) 
represent the number of distinct configurations of energy 
E in the j-th basin and introduce the entropy Sj{E) = 
In Wj {E). We now introduce the shifted PF 

z,{T)^ W,{E)e-^(^-''^ (7) 

EeAjE 

of the j-th basin and the free energy function 
fj{£j,E,T) = E - £.j - TSj{E), determined by 
the general summand in (0. The form of fj{£j,E,T) 
assumes that £j , E are independent, which is consistent 
with what was said above about the energy E of equilib- 
rium configurations and its independence from the basins 



to which they belong. The conditions for the minimum of 

fj {£j ,E,T)a.tE = 'E.J = E^ (T) are 

{dfj{£j,E,T)/dE),^^T\E^-^^=Q. (8a) 

{d^fj{£j,E,T)/dE^),^,T\^^-^^ > 0, (8b) 

and simplify to 

idSjiE)/dE),^\^^^^=P, dEj{T)/dT>0. (9) 

Both conditions are always met. The average basin energy 
Ej(T) determines the average entropy Sj(T) = Sj{E = 
Ej), so that fj{T) = -T In Zj{T) = Ej{T) - - 
TSj (T) represents the basin free energy. We wish to em- 
phasize that Ej{T), Sj{T), and fj{T) do not represent 
equilibrium quantities yet; the latter are determined only 
after Z{T) is evaluated. (If each basin is treated as repre- 
senting an independent system in a formal sense, then these 
quantities do represent equilibrium values for the particular 
basin.) 

We now group basins, indexed by j(A), into basin 
classes (BC) B\, indexed by A, so that all basins in a class 
have the same BM energy £ = £\. The basins in a class do 
not have to be close in the configuration space. The num- 
ber of basins in Bx is Nbm{£\), and the corresponding BM 
entropy is Sbm{£x) = ^nNBuiSx)- Let 

ZxiT) = ^liT), zx = Zx{T)/N^m{Ex). (10) 

denote the shifted and the mean shifted basin ;Sa-PF, re- 
spectively, so that 

Z{T) = ^e-^^^+^™(^^)zA(T), (11a) 

A 

£{T) = YS^e-^^"^'^'"''^'^Zx{T)/Z{T). (lib) 

A 

Here, £ = £(T) represents the equilibrium BM energy. It 
is easy to see that d£ / dT is a cross-correlation so that it 
need not have a unique sign 1 10]. The equilibrium free 
energy, entropy and energy are F{T) = —TlnZ{T), 
S{T) = -dF/dT and E{T) = F{T) + TS{T), 
respectively. 

SW Assumption. In the SW analysis, the sum over A 
in jllalllbll is replaced by a sum over the BM energy £x 
by assuming that Zx depends explicitly on £x in addition 
to T : Zx = Z\(£\ , T) . This issue has been examined 
earlier by us lioll . Here, we pursue the consequence 
of this assumption for the stability of this approach. Let 
fx{£x,T) = —Tlnzx{£x,T) be the mean free energy 
resulting from the mean basin PF Zx and Sx {£x , T) = 
— [dfx{£x,T)/dT]g^ the mean basin entropy. From the 
form of Zx (T) and Zx {£x ,T), it is obvious that we can 
rewrite Zx{£x,T) as follows: 

zx{£x,T)= Y VFA(^)e-''(^-^^)-^™(^^^ (12) 
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where Wx (E) represents the number of configurations of 
energy E that belong to Bx and indicated by £^ G AxE. 
Introducing Sx{£x,E) = lnWx{E) — SBui^x), we find 
that the general summand in J12ll determines a mean free 
energy function fx{£x, E,T) = E - £x - TSxi^x, E) 
whose minimization with respect to E at fixed £x , and T 
determines the mean free energy fx {£x , T) . Let Ex = 
Ex (T) denote the location of the minimum, the conditions 
for which are exactly the same as for fj{£j,E, T) above, 
except that the index j is replaced by A and Sx {£x , -E) is a 
two-variable function. Again, £x, and E are treated as two 
independent variables for the minimization to be carried 
out. The conditions for minimization are 

{dSx{£x,E)/dE)sAE^-E^ = /3, dEx{T)/dT > 0. 

(13) 

The mean free energy and entropy are given by f{£x,T) = 
fxi£x,Ex,T) = Ex - £x - TSx{£x,Ex) and 
Sx{£x,T) = S'a(£'a, -Ea)i respectively. It is easy to see 
that 

Sx{£x,Ex) = -[dfx{£x,T)/dT]^^ (14) 

as expected, where we must use 

{dEx/dT) =T[dSxi£x,Ex)/dT]e,, (15) 

which follows immediately from the first condition in (I13t . 
Since Ex is independent of £x at fixed T, we can differen- 
tiate fx{£x,Ex, T) with respect to £x to obtain 

{dS{£, T)/d£)T = + {df{£, T)/d£)^], (16) 

where fixed T means keeping Ex and T fixed simultane- 
ously, and where we have suppresed A to treat £ a variable. 

Zero-Slope Condition. Because of the assumed in- 
dependence of Zx{T), the general summand in ()llallllb|) 
becomes an explicit function of £, and we can mini- 
mize the corresponding free energy function Fb{£, T) = 
£ + f{£,T) — TSbm{£) with respect to £ at fixed 
T to determine Z{T). The minimum of Fb{£,T) is 
given by the conditions [dFB{£,T)/d£]T = 0, and 
[d^FB{£, T) /d£^]T > 0. The first condition is satisfied at 
the equilibrium BM-energy £= £{T) = £{T), see (lllbt . 
It is also given by the solution of 

dSBM{£)/d£ = P[l + {df{£, T)/d£)^], (17) 

and determines the equilibrium free energy F{T) = 
Fb{£,T), BM-entropy S'bm(7')_ = 'S'bm(^' = £), mean 
basin free energy fb(T) = f(£,T) and mean basin en- 
tropy 5b(r) = S(£,T) = -[df(£,T)/dT],; see (d. 
The equilibrium mean basin energy Eb (T) is obtained by 
the fundamental relation Eb{T) — £ = fh{T) + TS\,{T). 
It is easy to see that the form of the equilibrium free energy 
F{T) = f(£, T) +£{T) - TSbm(£) is the same as the 
free energy obtained by Goldstein in ©, except that the 
equations determining the equilibrium BM-energy are dif- 
ferent; compare © and l ll7ll . The two conditions become 



identical if / is taken to be independent of as was as- 
sumed by Goldstein. Comparing illi with jl6ll applied at 
£ = £,we obtain an interesting relation 

dSBM(£)/d£ + {dS(£, T)/d£)T = 0, (18) 

which will play a very important role in the following when 
we investigate the stability of this approach. 

The entropy S{T) can now be obtained by using the 
relation S{T) = -dF{T)dT. We immediately find that 
S{T) = S{£, T) + Sbm[£)- The equilibrium energy given 
by F{T)+TS{T), thus, turns out to be E^iT) introduced 
above. From the conventional analysis, this energy was 
identified as E{T). Thus, 

Eb{T) = E{T). 

Now, we apply O at^A = Eb{T) = E{T) to find 

{dE{T)/dT) =T[dS{£,E)/dT]j. (19) 

It should be noted that the entropy derivative on the right- 
hand side is the intrabasin change in the basin entropy with 
T without leaving the basin (fixed £). Comparing this with 
we find that the right hand side in both equations must 
be the same. This can only happen if jlSI l is fulfilled; we 
assume that [d£ /dT) ^ 0. This provides another justifi- 
cation for the validity of (I18t . and is merely a consequence 
of the first condition of stability. 

Curvature Condition. We now proceed to discuss the 
second condition for minimization. This condition of sta- 
bility dA£ = £ reads 

{d^f{£, T)/d£\ - T{d^SBM{£)/dt) > 0. (20) 

We differentiate (I17t at arbitrary £, which yields 

d'SBMi£) _ dl dfi£,T) 
d£^ d£^^^ d£ '^^^ 

d'f{£,T)dT^ 



dTd£ d£' 



We now set £" = f and use it in (I20t to finally obtain the 
condition to be 

[dSBM(£)/d£-d'f(£,T)/d£dT]{dT/d£)>0, (21) 

where we have used (I17t . Applying (1141 at equilib- 
rium, we obtain S(£,T) = -[df(£,T)/dT]j. Thus, 
the numerator in (I21t can be reduced to dSBM{£) / d£ + 
[dS{£, T) /d£]T- [The numerator can also be alternatively 
expressed as dS{T)/d£ - [dS(£,T)/dT]j/{d£/dT) 
where S{T) = S(£,T) + 5bm(^)]- Thus, the second 
condition of stability reads 

[dSBM(£)/d£ + [dS(£,T)/d£]T]{dT/d£) > 0, (22) 

which can never be satisfied in view of jlSIl unless 
{d£/dT) = 0. Since it is evident from (fTTEt that 
{d£ I dT) 7^ in general, we conclude that the curvature 
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of the free energy function at £ = £ must vanish on ac- 
count of the first condition of stability. Thus, we have fi- 
nally shown that the SW analysis is internally inconsistent 
and fails to provide a stable description of the system. 

Complexity Approach. We provide an alternative ap- 
proach ifioll which, as we show below, turns out to be 
a consistent and stable approach. We consider the un- 
shifted basin partition function Zj{T) = e^^^Zj{T) 
and introduce the unshifted basin free energy iPj{T) = 
— T In Zj (T) . The conditions of stability for the basin free 
energy ipj{E,T) = E — TSj{E) are given in (|9j. The 
basin free energy (pj (T) varies from basin to basin and rep- 
resents a family of functions, one for each j. Let J\f{ip, T) 
represent the number of basins having the same free energy 
93 for a given T and rewrite (0 as 

Z{T) = Y,M{^,T)e-f^. (23) 

The complexity is defined by S{ip, T) = In A/'((/?, T), in 
terms of which the conditions of stability at (p = Tp^^ = 

{dS{^,T)/d^)T\^^- = P, [d'S{lp„T)/d^l]T < 0. 

_ (24) 

The equilibrium complexity S{T) is given by S{^p^„T) 
evaluated at 93 = P\,- We consider the case so that S{ip,T) 
can be inverted at fixed T to express p as a function of 
S,T : if = ip{S, T). Thus, 

difiS, T) = (dif/dS)^ dS+ {dp/dT)^ dT. (25) 

At equilibrium, Tp^ = ip{S, T), and the coefficient of the 
first term becomes T according to the first relation in (1241 . 
In general, the coefficient of the second term is the negative 
basin entropy: S{ip,T) = —{d(p{S,T)/dT)s; compare 
with (I14t . Let us introduce cr{ip, T) = {dS{ip, T)/d(p)rp 
so that cr(^j,, T) = /? at equilibrium. From (I25t . we find 

{dSi^,T)/dT)^ = -dS{^,T)/dp)T {dp{S,T)/dT), 
= a{ip,T)S{ip,T). (26) 

We differentiate a{(p,T) with respect to T at constant (p 
and use the above equation to obtain 

{da{y^,T)/dT)^ = d^S{p,T)/dipdT 

= {d[a{^,T)S{^,T)]/d^)T, 

which is used to calculate 

da{^,T)/dT = {dai^,T)/d^)T[{d^/dT) + S{^,T)] 
+ a{p,T){dS{p,T)/d^)T. 

We now differentiate the first condition in (I24ll and use the 
above equation at equilibrium to obtain 

idai^^,T)/dlp^)T[{d^JdT) + S{^^,T)] 
= -P'-P{dS(Tp„T)/dlp^)T. 



The basin entropy S{(p, T) can also be expressed as 
^(iS, T) , so that the basin free energy function can be writ- 
ten as (^(5, T) = £^(T)—TS'(5, T). From this, we obtain 

{dcp{S,T)/dS)^ = -T{dS{S,T)/dS)^. 

At equilibrium (S =S), the left-hand side is equal to T 
from the first condition in i2Ai . Thus, at equilibrium, 

{dS(S,T)/dS)^ = -1. 

Since_ _{dS{^^„T)/dlp^)T 
{dS{S,T)/dSUdS/d^,) = -{dS{Tp^,T)/dRp^) = 
-[/3 + {dS/dT)-^{dT/dlp,)] = -f3[l + 

S^,{T){^T/^(p^^)], where we have we used ( I26I I at 
equilibrium [a{Tp^^,T) = (3, and Sb(T) = S{'ip^^,T)]. 
Thus, we find that 

{da{Tp,,T)/d^,)r[{d^JdT) + S,{T)] 

= P'S,{T){dT/dip,), 

so that the second condition of stability becomes 

ida{Tp„T)/dlp,)r = P'S,{T){dT/dip^)/ {dS/dT) 
<0, 

where we have used the relation 

T {dS/dT) = d^JdT + 5b (T) (27) 

obtained from Therefore, the second condition of sta- 
bility finally becomes 

{dS/dT) {d^JdT) < 0. 

We expect (^dS/dT) to be positive at low temperatures, 
so the stability condition there reduces to {dip^,/dT) < 0. 
Because of (I27t . it is easy to see that (dS/dT^ > even 
if {d^p^^/dT) changes sign, so the stability always requires 

idS/dT) > 0; {dpJdT) < 0. 

It is our pleasure to thank Fedor Semerianov for various 
discussions. 
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